Glass glass transition and new dynamical singularity 
points in an analytically solvable p-spin glass like model 
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We introduce and analytically study a generalized p-spin glass like model that captures some of 
the main features of attractive glasses, recently found by Mode Coupling investigations, such as 
a glass/glass transition line and dynamical singularity points characterized by a logarithmic time 
dependence of the relaxation. The model also displays features not predicted by the Mode Coupling 
scenario that could further describe the attractive glasses behavior, such as aging effects with new 
dynamical singularity points ruled by logarithmic laws or the presence of a glass spinodal line. 

PACS numbers: 64.70.Pf, 05.50.+q 64.60.Ht 
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Recent Mode Coupling Theory (MCT) 1] investiga- 
tions discovered 2] a new kind of glasses, the so called 
attractive glasses, characterized by the presence of a 
glass/glass transition line where higher order dynamical 
singularity points are located. As MCT calculations trig- 
gered intense experimental researches 0, Q , their limit is 
that, by definition, they can correctly describe only the 
region outside the glass. In this paper, we introduce a 
schematic p-spin glass like model exhibiting a glass/glass 
transition which we can analytically investigated to es- 
tablish the proper scenario of the glassy phase. This leads 
to new testable predictions on the aging dynamics in that 
region, where new types of singularities are found. 

In MCT attractive glasses the potential between the 
particles presents attraction for some range besides the 
hard core repulsion Its volume fraction - temperature 
phase diagram shows two branches of the glass transition 
line, one due to the repulsive part of the potential, ex- 
tending to high temperatures for volume fractions near 
the close packing value, and one related to the attrac- 
tion, extending to low volume fractions up to meet the 
liquid/gas spinodal. For sufficiently narrow widths of the 
potential well one of the two branches interestingly enters 
in the glass region giving rise to a glass/glass transition 
line, ending in a higher order dynamical singularity point. 
This special point in the MCT framework is character- 
ized by a logarithmic time dependence of the relaxation. 
The theory applies to attractive colloids; actually, some 
evidence of the MCT results has been found by recent ex- 
periments U 0. Nevertheless, for some aspects, the ap- 
proach used in [2| turns out to be inadequate. It is known 
in fact that MCT is valid only for equilibrium dynamics 
and cannot be trivially extended inside the glassy region, 
where the glass/glass line itself is located and aging phe- 
nomena are present. Another problem arises when the 
phase coexistence interferes with the glass transition line. 
In this case MCT does not allow a self-consistent de- 
scription of the structural relaxation and the gas/liquid 
spinodal line is usually obtained a priori, not taking into 
account the presence of the glassy phase. The purpose of 



the present work is to begin to cure these inconsistencies. 
Exploiting the strong analogies between glassy systems 
and the class of 'discontinuous' spin glasses noticed in the 
last years p, we introduce a schematic solvable model 
providing a coherent picture of the glass/glass line, its 
off-equilibrium dynamics and singularity points. 

The model consists of a combination of two (or more) 
»-spin models 0, 0, E| diluted through density variables 
|3 • We analytically solve its Langevin dynamics and 
for comparison with MCT inside the glassy region we 
also use the Mode Coupling approximation. We show 
that: (A) The glass/glass transition line does not coin- 
cide with that predicted by MCT but follows a different 
curve (Figs. 1,2). (B) Along the glass/glass line (see 
Fig. 1) in the early regime the stationary autocorrela- 
tion function $(t — t') exhibits, like in MCT, a power law 
relaxation except at its endpoint As where is logarith- 
mic. (C) In the long time regime, where MCT predicts a 
plateau, we find an interesting aging dynamics. In par- 
ticular, two more dynamical singularity points, B 3 and 
B' 3 (one for each kind of glass, see Fig.l) characterizes 
the relaxation in the aging regime. On one side of the 
glass/glass line from C to B 3 (resp. B' 3 on the other side, 
see inset in Fig. 1) the long times relaxation is a power 
law, except right at B 3 where it is logarithmic (analo- 
gously for B' 3 ). From B 3 (resp. B 3 ) to A 3 , the aging 
behavior is spin glass-like 5|. (D) Finally, at the A^ sin- 
gularity (i.e., in the case where ^3 coincides with point 
C in Fig. 1) we show that the logarithmic decay of the 
correlation function in the stationary regime $ (t — t') is 
followed by an aging regime ruled again by a logarithmic 
law. (E) We calculate in a consistent way the density- 
temperature phase diagram and the spinodal lines (see 
Fig. 1) and find they are influenced by the presence of 
the glassy phase. 

The Model. We consider a lattice gas model, where 
particles (rij = 0, 1) interact in groups of p via their spins 
(sj = ±1) through quenched random couplings J^...^, 
with zero mean and variance p\Jp/ (2iV p_1 ), and also in 
groups of r through (nonrandom) coupling constants K r 
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(typically below we take two or three values of p and one 
value of r). The mean field Hamiltonian is given by: 
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where is the chemical potential. The spherical version 
of the model Q is here considered: we introduce two 
new spin variables s a i (a = 1,2), defined by Si = s%i, 
n-i = {suS2i + 1) /2, and enforce for these ones spheri- 
cal constraints s 2 ^ = TV) [T^j. To exclude the pos- 
sibility of continuous transitions, not related to a glass 
transition, the restriction p > 2 is taken |(J. The two 
spin fields s a i (a = 1,2) evolve via usual Langevin equa- 
tions with thermal noises £ a i (t) , having zero mean and 
variance (£ a , (t)£bj (t')) = 2TS a bSijS (t — t'). Due to its 
quadratic nature the model is exactly solvable. It has 
a gas, a liquid and a glassy phases, respectively called 
P~ , P + and G (see Fig. 1). We focus here on the 
dynamics, and in particular on the correlation func- 
tion C (t, if) = J2 ia ( s u (<) s <« (*')) / (2N), the related re- 
sponse function G{t,t') and the density d(t) — C (t,t). 
Standard functional techniques allow to derive the 

following dynamical equations for these quantities after 
a rapid quench at t — from high temperature (t > t'): 



dC (t,t') 
dt 



z(t) 



C(t,t') + 2TG (t',t) 



+- / dutp' (C(t,u))G(t',u) 



+- du<p"(C(t,u))G(t,u)C(u,t'p) 



dG (t,f) 
dt 



+ - / duip" (C(t,u))G(t,u)G(u,t'X3) 
2 Jt> 

d' (t) = (z (t) + 2^ + 2 X ' (d (t))) (1 - d (t)) - T (4) 

Here we define <p (x) = (1/2) JpX p , x( x ) — 
(l/2)Er K rX r and z(t) = Zl (t)-n - x' (d(t)) with 
Z\ (t) = Z2 (t) the two time-dependent Lagrange mul- 
tipliers related to the spherical constraints (the prime 
stands for the first derivative) . The above equations com- 
pletely define the dynamical model, given the initial den- 
sity d(0). In the following we consider the large times 
limit, where one-time quantities reach stationary values: 
d = d(t — > oo ) and z = z (t —> oo). 

MCT calculations. Assuming equilibrium dynamics, 
we have that C = C(i-t') and G = G{t-t') Vt, t! and the 
Fluctuation-Dissipation theorem (FDT) holds: TG (r) = 
— C (t) (with r = t — t'). Under these hypotheses, Eqs. 




FIG. 1: The model dynamical phase diagram in the density- 
temperature plane (the parameters are here: pi = 3, p2 = 11, 
j| = .74, jf t = 4 and r = 12, = —12). One can have a 

fluid (P) phase consisting of gas (P _ ) and liquid (P -1- ), or a glass 
(G) phase. Continuous lines depict the liquid/glass (P + /G) and 
glass/glass (G/G) transition lines, the spinodal of the fluid (P~-sp 
and P+-sp) and that of the glass (G-sp). The inset shows an en- 
largement of the region around the G/G transition line. The dotted 
line represents the incorrect G/G line predicted by MCT. On the 
G/G line (full line from C to A3), along with its endpoint A3, there 
are new dynamical singularity points, P3.P3, characterized by a 
logarithmic relaxation. Also the crossing point G between the two 
branches of the P+/G transition line is shown. 

(L'IMI) yield the long time schematic MCT equation ^1 
for the correlator $ (r) = C (r) /d, 



To—r^ + $ (t) + 

dr 



dT'm{T-T , )^j^l = (5) 
dr 



where m (r) = <p' (C (r)) d/T 2 . In the standard MCT no- 
tation |l| the kernel m (t) is usually written as m (r) = 
F ($ (t) , v) with F (/, v) = J2 n v nf n -, this allows to es- 
tablish a mapping relating T, d, J p and the usual MCT 
control parameters v n , given by v n — pJ 2 d p / (2T 2 ) 
where n = p — 1. If p\ 1 p2 1 . ■ ■ denote the values as- 
sumed by p, we deal with the so-called F Pl -i,p 2 -i.... MCT 
model. The non-ergodicity parameter / = $ (r — > 00) is 
obtained by solving the bifurcation equation Z (/) = 1 
(i.e., the r — > 00 limit of Eq. JSJ), where 



£(/) = 



1 



1-/ 



(/,<;) = 



1 



1-/ 



2T 2 



(6) 



The inequality 1 < Z ($ (r)), deriving from the decreas- 
ing character of 3>(t), implies that f is the largest so- 
lution of the bifurcation equation [lj. The conditions 
Z (f c ) = 1 and Z' (f c ) = determine the transition line 
in the plane d — T (see Fig. 1). The P + /G transition, 
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where / jumps from zero to a nonzero value, f c , corre- 
sponds to singularity points of type A 2 Q. MCT Eq. 
(J5J describes the relaxation when this critical line is ap- 
proached from the P + phase: $(r) decays in two steps 
characterized by power laws with critical exponents de- 
termined by the parameter A = 1 — (1 — f/f Z" (f c ) /2 < 
1 0. Inside the glassy region, for a set of values of the 
parameters v n , an other line of A 2 singularities appear 
(dotted line in the inset of Fig. 1) where two non zero 
solutions for the dynamical order parameter, /, are al- 
lowed. This line defines the glass-glass {G/G) transition, 
where a first glass./vi , transforms discontinuously into 
a second one, f C 2 [H El H^. Its endpoint, where A = 1, 
i.e., Z" (f c ) — 0, results to be an higher order singularity 
A3 (with f c i — f c2 ). To have a G/G line in our model, 
at least two values of p are needed Q: as an example, 
for pi — 3 the lowest other possible value is p 2 = 11 (the 
case of Fig.s 1,2). 

The glassy region and the new G/G line. Eq.J^J), i.e., 
MCT, is correct only out of the glassy phase (for / = 0), 
where the equilibrium dynamics assumption holds. Thus, 
the above results on the G/G line are not correctly de- 
scribed by this theory. Within our model, we discuss now 
the more complex non-equilibrium equations governing 
the glassy phase, the correct position of the G/G transi- 
tion and the new aging behavior along the G/G line. In 
the glassy phase, as in the usual p-spin glass model, the 
relaxation of <&(t, t') can be split in two parts 
in the first, describing the so called FDT regime (t ~ t'), 
<f>(t,t') is a function, $fdt(t), of the times difference 
t = t — t'; the second part for t 3> t' corresponds to an 
aging regime described by a function ^^(i, t'). In the 
FDT regime, i.e., the approach of the correlator to the 
plateau f c , $fdt{t) satisfies the following equation: 



to - 



d$ (r) 
dr 



+Z$ (t)+(1 - Z)+ [ dr' 
Jo 



w(T " r) ^^ =0 
(7) 

where the existence of a well defined aging solution fixes 
Z to the minimum value (not necessarily equal to 1) sat- 
isfying the stability requirement Z (f) < Z (&fdt (t)), 
i.e., the one corresponding to the marginal stability con- 
dition Z' (/) = @, [H EH (see Fig. 2). Note that 
Eq.0 coincides with © only for Z = 1. The G/G 
transition line results thus to be modified with respect 
to MCT and located by the conditions Z (f c i) — Z (f C 2) 
and Z' (f c i) = Z' (/ c2 ) = (see Fig.s 1,2). Its endpoint 
A 3 (see Fig. 1), is characterized, like in MCT, by one 
more condition: f c \ = f c2 = f c with Z" (f c ) — 0, corre- 
sponding to A = 1. 

In the aging region, the correlator Cag (*>*') = 
^ ag (t, t') d and the response function Gag (t, t') 
obey a generalized FDT relation, TGAG(t,t') — 
xOCag (t>t') /dt', where x < 1 is a constant with the 
physical meaning of the ratio between the bath temper- 
ature and an effective temperature T e tf, describing the 
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FIG. 2: The function Z (/) denned in Eq. J|J is plotted for several 
decreasing temperatures T (from the top to the bottom) and fixed 
density. For high T the liquid solution (P + ) is given by the value 
/ = corresponding to the lowest value of Z namely Z = 1. By 
lowering T the glass transition T c is reached when the minimum of 
Z(f) approaches the level Z = 1 (bold line marked P+/G). For 
T < T c the plateau f c is given by the value of / where Z(f) has 
the minimum. For low enough T, Z(f) develops a second mini- 
mum. The G/G transition is reached when the two minima have 
the same depth (bold line marked G/G). Note that in MCT ap- 
proximation this would be reached instead when the second mini- 
mum approaches the level Z = 1 (line marked G /G MCT ). 

out of equilibrium configurations x = T/T e ff 0. 

The dynamics along the G/G line. On the G/G line 
two glasses coexist corresponding to two different val- 
ues of the plateau, /, and thus two different values of 
x = x(f) and A = A(/) (with < x < 1 and < A < 1), 
except at the endpoint A3 where the glasses coincide. 
As shown below, along the line the value of A deter- 
mines the properties of the dynamics in the FDT regime, 
namely the approach to the plateau, while the ratio A/x 
determines the aging regime, i.e., the departure from the 
plateau. Four different behaviors are found along the 
G/G line (see Fig. 1): the case X/x < 1, corresponding 
to the points from C to B3 (resp. B' 3 , on the other side of 
the line, see below); the new dynamical singularity points 
B 3 and B' 3 where X/x = 1; the case X/x > 1 extending 
from S3 (resp. B' 3 ) to the endpoint A3; and the singular- 
point A3 itself where X/x > 1, but A = 1. 

More precisely the approach to the plateau in the 
FDT regime, along the entire G/G line, is given by 
<&fdt (t) — f c ~ , where the exp onent j3 is given by 
T 2 (1 - /3) /r (1 - 2/3) = A d S Ef (r is Euler gamma 
function). At the endpoint A3 where A = 1 this is re- 
placed by a logarithmic behavior given by <&fdt (t) — 
/ c ~l/ln 2 r 

The departure from the plateau in the aging regime, 
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along the G/G line from the crossing point C up to the 
two new points B3 and B 3 , one for each side of the line, 
is given by a power law: 

<S>A G (t' +r,t')- f c ~-(r/T(t>)) a (8) 

with r 2 (l + a)/r(l + 2a) = X/x &}. Here the t'- 
dependence is contained in T (t r ) ~ if 

At the two new dynamical singular points, B3 and B 3 , 
where X/x = 1, this power law is replaced by a logarith- 
mic behavior, one for each side of the G/G line. For in- 
stance, by approaching B 3 from the right (the side where 
the plateau value f c is higher) one has 

$A G (t' + r,t')-/c~l/ln 2 (r/T(t')) . (9) 

From the point B3 to the endpoint A3 (analogously 
from B' z to A3 on the other side of the G/G line), 
X/x > 1. This implies that the present one-step Replica 
Symmetry Breaking solution should not hold an d, 
instead, a spin glass-like aging dynamics should be found 
in such a region 

Interestingly the length of the G/G line can be varied 
and, in a version of the model with three p- values [lij . 
the A3 endpoint singularity can be made to coincide with 
the crossing point C (see Fig. 1). In this particular case 
the A3 becomes an A4 singularity, as also found by MCT 
,2\ (the simplest model with such a singularity is found 
to be pi = 3,p 2 = 4,p 3 = 11 14]). In this case the FDT 
regime is logarithmic $fdt (t) — fc ~ 1/ l nr , a s in MCT 
ETUI EH > as we H as the aging regime: 

$ AG (f + r, - fc ~ 1/ In (r/T (f)) ■ (10) 

Spinodal lines. We show now that the glass transition 
curve P + /G intersects the fluid spinodal line determining 
the existence of a glass spinodal line. The two spinodals 
are obtained by enforcing the vanishing local stability 
conditions in Eq. (@J: 

/?y (d) + 2I3 X " (d) y 2 = (11) 

(l-d) 2 d 2 {l-f) 2 

The P~/P+-spinodal corresponds to a solution, T(c£), 
with / = 0, as for the G-spinodal solution the value of 
/ of the glassy phase is given by the marginal stabil- 
ity condition Z' (/) = 0. The two lines, completing the 
dynamical phase diagram, together with their range of 
validity are shown in Fig. 1. 

In conclusion, the model we have introduced is a sim- 
ple extension of p-spin glass models 0, f° r the glass 
transition. Our model exhibits a glass/glass transition 
line which can be analytically investigated. Actually, it 



turns out to be different from the G/G line derived by 
equilibrium MCT and to be characterized by new dy- 
namical singularity points and off-equilibrium properties. 
These results can help to shed deeper light on the prop- 
erties of attractive glasses, such as colloidal suspensions. 
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CRdC-AMRA, INFM-PCI, EU MRTN-CT-2003-504712. 
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